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INTRODUCTION 


With the advent of finite element models of varying complexity the focus 
of solutions to problems in solid mechanics has shifted very strongly into the 
direction of more accurate material description. This is especially true for 
materials for which strength characteristics vary widely with state of stress. 

In particular, concrete which is non-isotropic at any level of deformation and 
is also non-linear in terms of stress-strain relationships has been singled out 
for intensive study. This includes work on constitutive relations (refs. 1 to 
5), and failure (refs. 3,4, and 6 to 15). 

Related developments in failure theories are those included in references 
16 to 25. 

This list is by no means exhaustive but formed a background basis which led 
to the model proposed in this paper. In particular, the developments in the 
areas of maximum deformation theory and the Von Mises-Hencky theory provided 
motivation for the concept used here (ref. 26). 


FAILURE SURFACE 


The proposed generalized constitutive equation is an extension of the work 
of Hu and Swartz (26) on the study of the failure of materials so that for any 
kind of material in any state, of stress its mechanical behavior can be charac- 
terized by a single functional. According to the theory, a material failure 
initiates when the state of stress at a point is such that the following func- 
tional reaches a threshold value F : 
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In this equation 


a = a material related scalar factor which is determined by experi- 
mental data. 


J 1 = Cl + o 2 + o 3 . 

J 2 = (o-l - o 2 ) 2 + (o 2 - a 3 ) 2 + (a 3 - a^ 2 . 


°1» °2 ’ a 3 


the components of principal stresses, tensile stress has positive 
value. A vector of principal stresses is expressed by a. 
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the absolute value of the ultimate stresses of the material under 
uniaxial tension and compression, respectively. 
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t 
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the absolute values of strain component in the direction of uniaxial 
force when stress reaches the corresponding ultimate value. 


E(a) = the secant modulus of elasticity of the material subjected to uni- 
axial stress (tension or compression as appropriate). 

y = the Poisson's ratio of the material (related to stress level). 


CONSTITUTIVE EQUATIONS 


Non-linear response will take place when the increment of the functional is 
in an increasing manner and its value is beyond some threshold level (possibly 
zero) . The material will be fractured at points where the state of stress 
reaches the surface of failure, i.e., F(a) = 1. 

Between the initial state to fracture, the response is assumed to be charac- 
terized by Drucker's (17) theorem of orthogonality with the use of the functional 
proposed by the generalized failure theory, eqn. (1). That is, if the increment 
of principal strain components is decomposed into linear and non-linear 
increments , 

{de} = {de e } + {de P }, (2) 

the non-linear part is characterized by 

{de P } = G(a){g(o)} dF. (3) 


where 


G(a) = a scalar function of stresses, 

{g(o)} = a unit vector of the gradient of the functional at the 
point of interest. 

Note that G values (under uniaxial loading) can be determined by some stress- 
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strain relationship. In the general case, especially for those materials like 
glass or concrete, the G function in uniaxial tension differs from that of uni- 
axial compression remarkably. Therefore a weighted average is proposed. Accord- 
ing to this, the relationship between the vectors and principal strain increment 
and the principal stress increment is proposed to be characterized by the follow- 
ing generalized constitutive equation: 




l -y -y 


{de} = 

1 

E 

-y l -y 

+ G(ct) (g(a) } (_VF (a) J 



-p -y l , 



(4) 


In this E = the modulus of elasticity of the material, ^_VF | = the row vector of 

the gradient of the functional at the point of interest, and the function G(ct) 
is calculated by 


G(a) 



H(F-F ) H(AF) 
o 


(5) 


which is a weighted average according to the values of the corresponding princi- 
pal stress components. G (cr^) depends on the value of the i-th uniaxial princi- 
pal stress and is calculated according to 


G i (a i } = G+ f F(CT i>J H(a i) + G EF(ct ± )] H(-a ± ). 


( 6 ) 


In this formula, H is a unit step function. 


AF = F(a. 
Then, 1 
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(7a) 


(7b) 


where g^ is the i-th component of the unit vector {g}. 

C £ t Q 

If a ^ = a f » = e f and “ = 1 is selected, the proposed failure theory 

agrees with the Von Mises ' theory and the generalized constitutive equation 
reduces to the well known Prandtl-Reuss stress-strain relation (refs. 4,11). 
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NUMERICAL EXAMPLE 


Applying the proposed theory to plain concrete for the purpose of illustra- 
tion equ. (1) is modified to read 


F 


a_ 
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c t , c t. , 
°£ °f - (<, f -°£ )J ] 
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[o^y (cr 2 +cr 3 ) ] (o f t +o f c ) 
a f t a f c (l+y) + (a f t - V o f C )J 1 


( 8 ) 


The implication used is an invariant modulus of elasticity and for concrete, so- 

called initial yield occurs at about .45 f' . 

c 

Using test data for concrete (3) the level surfaces of the functional, the 
stress-strain curves of uniaxial tension and compression tests and the variations 

of G + [F(a^)] and G [F(o^)] are shown in Figs. 1-3. The best value of a for 

concrete is 0.46 which has been used in these curves. 


Note in Fig. 1 the shape of the curve in the tension-compression zone follows 
very closely the shape of the experimental curve obtained by Kupfer, Hilsdorf 
and Rusch (3). The stress-strain curves presented in Figs. 2 and 3 were obtained 
using testing equipment described in Ref. 27. Using these data, Equations 8 for 

F and 7a and 7b for G + and G~ were evaluated numerically to obtain the curves 
displayed. 


CONCLUSIONS 


A proposed constitutive model for non-linear materials has been presented. 

The primary virtues of the model are its logical combination of distortion states, 
inherent simplicity and generality. 

Results presented for the model applied to concrete show good agreement with 
published experimental data for failure. The model can be readily incorporated 
into existing computer codes provided sufficient experimental supportive data 
are available. 

As reported elsewhere (10), for instance, the experimentally obtained para- 
meter a for cast iron is around 0.93. 

The proposed constitive equation is presently being utilized in the develop- 
ment of a finite element code for determination of unstable crack growth and 
stress intensity in concrete beams. 
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Figure 3.- Data of uniaxial compressive test of concrete. 
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